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ABSTRACT 

In this paper, I revisit the constraints obtained by several authors (Reichart 
et al. 1999; Eke et al. 1998; Henry 2000) on the estimated values of fi^, n 
and (Ts in the light of recent theoretical developments: 1) new theoretical mass 
functions (Sheth & Tormen 1999, Sheth, Mo & Tormen 1999, Del Popolo 2002b); 
2) a more accurate mass-temperature relation, also determined for arbitrary Vt^ 
and (Voit 2000, Pierpaoli et al. 2001, Del Popolo 2002a). Firstly, using the 
quoted improvements, I re-derive an expression for the X-ray Luminosity Func- 
tion (XLF), similarly to Reichart et al. (1999), and then I get some constraints 
to fim and ra, by using the ROSAT BCS and EMSS samples and maximum- 
likelihood analysis. Then I re-derive the X-ray Temperature Function (XTF), 
similarly to Henry (2000) and Eke et al. (1999), re-obtaining the constraints on 
n, CTg. Both in the case of the XLF and XTF, the changes in the mass func- 
tion and M-T relation produces an increase in Vt^ of ~ 20% and similar results 
in (Tg and n. 

Subject headings: cosmology: theory - large scale structure of universe - galaxies: 
formation 



1. Introduction 

Galaxy clusters represents the virialization stage of exceptionally high peaks of initial 
density perturbations on comoving scales of ^ lO/i^^Mpc, and as such they provide useful 
cosmological probes. The evolution in the abundance of clusters is strongly dependent on 
the cosmological density parameter, (Evrard 1989; Oukbir & Blanchard 1992; Eke, Cole 
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& Prenk 1996; Donahue et al. 1998; Borgani et al. 1999): in a critical universe with 
Q = 1 , perturbation growth continues forever, while in a low-density universe (fl < 1), 
growth significantly decelerates once z ~ ^2^^ — 1. Although the quoted cosmological test 
is in principle very powerful, there are two main problems in practical applications. First, 
theoretical predictions provide the number density of clusters of a given mass, while the 
mass itself is never the directly observed quantity. Secondly, it is needed a cluster sample 
that spans a large z baseline, and is based on model- independent selection criteria ^. In this 
respect. X-ray observations provide a very efficient method to identify distant clusters down 
to a given X-ray flux limit, and hence within a known survey volume for each luminosity, 
Lx- For this reason, most studies using clusters as cosmological probes are based on X-ray 
samples. It is well known that clusters are strong X-ray emitters and so cluster evolution 
can be inferred from the study of X-ray properties of distant clusters. The amount of 
observational data concerning high-redshift cluster properties has increased in the past years. 
EMSS (Henry et al. 1992; Gioia & Luppino 1994), ASCA (Donahue 1996; Henry 1997) and 
ROSAT (Ebeling ct al. 1997; Rosati ct al. 1998) measurements of the X-ray emitting 
intracluster plasma complements low-redshift studies carried out by Edge et al. (1990) 
and David et al. (1993). In addition, galaxy velocity dispersions for a well-defined sample 
of high-redshift clusters (Carlberg et al. 1996) are provided by the CNOC survey. The 
Press & Schechter (1974) (hereafter PS) formalism has been heavily used to model the 
cluster population. The combination of the PS mass function and the X-ray cluster catalogs 
represents a unique opportunity to constraint cosmological parameters, (e.g. the mass density 
parameter, Om)- The PS approach has gained with time an increasing favor since this 
approach ofi^ers a number of advantages when compared with more traditional methods of 
measuring fl,^. Firstly, unlike methods that only probe Qm over small spatial scales, the PS 
approach probes over the scales of the X-ray cluster catalogs. Secondly the PS approach 
seems to be relatively insensitive to the presence of a cosmological constant (Henry 1997; Eke 
et al. 1998; Viana & Liddle 1999). Thirdly, till some years ago, it was shown that numerical 
simulations reproduce the PS mass function quite accurately (Eke et al. 1996; Bryan & 
Norman 1998). More recent studies has, however, shown some discrepancies between PS 
and simulations. Although the analytical framework of the PS model has been greatly 
refined and extended (e.g., Lacey & Cole 1993), it is well known that the PS mass function, 
while qualitatively correct, disagrees with the results of N-body simulations. In particular, 
the PS formula overestimates the abundance of haloes near the characteristic mass M* and 
underestimates the abundance in the high mass tail (Efstathiou et al. 1988; Lacey & Cole 
1994; Tozzi & Governato 1998; Gross et al. 1998; Governato et al. 1999). The quoted 
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discrepancy is not surprising since the PS model, as any other analytical model, should 
make several assumptions to get simple analytical predictions. 

There are different methods to trace the evolution of the cluster number density: 

a) The X-ray temperature function (XTF) has been presented for local (e.g. Henry & Ar- 
naud 1991) and distant clusters (Eke et al. 1998; Henry 2000). The mild evolution of the 

XTF has been interpreted as a strong indication for a low density universe (0.2 < < 0.6). 
As described by Colafrancesco, Mazzotta & Vittorio (1997), Viana & Liddle 1999, and Blan- 
chard et al. 2000, this conclusion could be weakened by uncertainties related to the limited 
amount of high-z data and to the lack of a homogeneous sample selection for local and dis- 
tant clusters. 

b) The evolution of the X-ray luminosity function (XLF). In this case, we need a relation be- 
tween the observed Lx and the cluster virial mass. Studies of the intra-cluster-medium (ICM) 
(Ponman, Cannon & Navarro 1999) showed that non-gravitational heating and, possibly, ra- 
diative cooling significantly affects the Lx — M relation. At the same time, other studies 
(Borgani & Guzzo 2001) showed the X-ray luminosity to be a robust diagnostic of cluster 
masses. Furthermore, the most recent flux limited cluster samples contain a large (~ 100) 
number of objects, which are homogeneously identified over a broad resdshift baseline, out 
to z ~ 1.3. 

The results for obtained are several. 
Kytayama & Suto (1997) and Mathiesen & Evrard (1998) analyzed the number counts from 
different X-ray flux-limited cluster surveys and found that the resulting constraints on are 
rather sensitive to the evolution of the mass luminosity function. Analyzing EMSS, Sadat, 
Blanchard & Oukbir (1998) and Reichart et al. (1999) found results consistent with O,^ — 1. 
A result consistent with Qm ~ 1 was found by Blanchard & Bartlett (1998), and Viana & 
Liddle (1999) found that ^ 0.75 with > 0.3 at the 90% confidence level and ^ 1 
still viable. Blanchard, Bartlett & Sadat (1998) found almost identical results (Jim — 0.74 
with 0.3 < r^in < 1.2 at the 95% confidence level). Eke et al. (1998) found = 0.45 ± 0.2. 
It is interesting to note (as previously mentioned) that Viana & Liddle (1999) used the same 
data set as Eke et al. (1998) and showed that uncertainties both in fitting local data and in 
the theoretical modeling could significantly change the final results: they found ~ 0.75 
as a preferred value with a critical density model acceptable at < 90% c.l. 

Different results were obtained by Bahcall, Fan & Cen (1997) {Q^i = 0.3 ± 0.1), Fan, 
Bahcall & Cen (1997) (^n - 0.3 ±0.1), Bahcall & Fan (1998) {n^ = 0.2+°;?). The previous 
example together with other not cited, shows results span the entire range of acceptable 
solutions: 0.2 < < 1 (see Reichart et al. 1999). 

The reasons leading to the quoted discrepancies has been studied in several papers (Eke 
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et al. 1998; Reichart et al. 1999; Donahue & Voit 1999; Borgani et al. 2001). According to 
Reichart (1999) unknown systematic effects may be plaguing great part of tlie quoted results. 
Systematic effects entering the quoted analysis are: 1) The inadequate approximation given 
by the PS (e.g., Bryan & Norman 1997). 2) Inadequacy in the structure formation as 
described by the spherical model leading to changes in the threshold parameter 6c (e.g., 
Governato et al. 1998). 3) Inadequacy in the M-T relation obtained from the virial theorem 
(see Voit & Donahue 1998; Del Popolo 2002a). ^ 4) Effects of cowhng flows. 5) Determination 
of the X-ray cluster catalog's selection function. 

Voit & Donahue (1999) point the attention on similar and different items to that stressed 
by Reichart et al. (1999): 1) Deviation from the Press-Schechter orthodoxy (similarly to 
Reichart et al. 1999). 2) Missing high redshift clusters in the data used (e.g., the EMSS). 
3) Inadequacy in the M-T relation. 4) Evolution of the L-T relation. 

Eke et al. (1998), in a detailed study of the systematic uncertainties in the determination 
of Vtui showed that even the differences between the derived best-fitting parameters from the 
and the Maximum likelihood techniques has non-negligible contributions in the systematic 
uncertainties. As stressed by Henry (2000), the maximum hkehhood fit only provides the best 
fit, it does not provide an assessment of whether that fit is a good fit to the data. Marshall 
et al (1983) require that the two independent cumulants (sec Eqs. 7, 8 of Henry (2000)) 
both be uniform according to the two-tailed Kolmogorov-Smirnov test. They recommend 
rejecting the model if the probability of either being uniform is less than 0.05. 

Although the quoted uncertainties has been so far of minor importance with respect to 
the paucity of observational data, a breakthrough is needed in the quality of the theoretical 
framework if high-redshift clusters are to take part in the high-precision-era of observational 
cosmology. 

Even if there have been several recent and detailed studies of the cluster abundance 
(Bond & Myers (1996), Eke, Cole & Frenk (1996), Viana & Liddle (1996; 1999), Co- 
lafrancesco, Mazzotta & Vittorio (1997), Kitayama & Suto (1997), Eke et al. (1998), Pen 
(1998), Wang & Steinhardt (1998), Donahue & Voit (1999) and Henry (2000)) 

The above discussion and recent developments in terms of both theory and observa- 
tion which suggest it would be useful to revisit this constraint: 1) the addition of ASCA 
temperatures (Tanaka, Inoue & Holt 1994) means that there is now a well defined local 
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temperature function for clusters, with relatively small errors in temperature. 2) A num- 
ber of large N-body simulations have accurately determined the mass function of virialized 
haloes (e.g. Governato et al. 1999), finding non- negligible deviations from the old Press- 
Schechter (1974; hereafter PS) theory. 3) More ambitious hydrodynamical simulations of 
cluster formation (e.g. Frenk et al. 1999, and references therein) and theoretical analysis 
(Voit 2000, Del Popolo 2002a) have resulted in improvements in the relationship between 
mass and temperature and a better estimate of its scatter. 

These reasons lead mc to rc-calculate the constraints on Q„^, n and ag, using the XLF 
and XTF. The paper is organized as follows: in Sect. 2, 1 rc-calculatc the XLF, as done by 
Reichart et al. (1999) and obtained constraints for VL^ and n using data from ROSAT'BC'^ 
and EMSS samples. In Sect. 3, I re-calculate the XTF, as done by Henry (2000) and Eke 
et al. (1998) and obtained constraints for and n and a^. Sect. 4 is devoted to results 
and to conclusions. 



2. Constraints to Vt^ and n from the XLF 

In this section, similarly to Reichart et al. (1999), I derive an expression for the XLF 
(using now the mass function and M-T relation obtained in Del Popolo 2000a, Del Popolo 
2000b, respectively) and then I get some constraints to Q,^ and n, by using the ROSATBC^ 
and EMSS samples. 

As previously quoted, although the analytical framework of the PS model has been 
greatly refined and extended, it is well known that the PS mass function, while qualitatively 
correct, disagrees with the results of N-body simulations (Efstathiou et al. 1988; White, 
Efstathiou & Frenk 1993; Lacey & Cole 1994; Tozzi & Governato 1998; Gross et al. 1998; 
Governato et al. 1999). The quoted discrepancy is not surprising since the PS model, as any 
other analytical model, should make several assumptions to get simple analytical predictions. 
As previously reported, the main assumptions that the PS model combines are the simple 
physics of the spherical collapse model with the assumption that the initial fluctuations 
were Gaussian and small. On average, initially denser regions collapse before less dense 
ones, which means that, at any given epoch, there is a critical density, 5c{z), which must be 
exceeded if collapse is to occur. In the spherical collapse model, this critical density does not 
depend on the mass of the collapsed object. Taking account of the effects of asphericity and 
tidal interaction with neighbors, Del Popolo & Gambera (1998) and Sheth, Mo & Tormen 
(2001) (hereafter SMT), using a parametrization of the ellipsoidal collapse, showed that the 
threshold is mass dependent, and in particular that of the set of objects that collapse at 
the same time, the less massive ones must initially have been denser than the more massive. 
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since the less massive ones would have had to hold themselves together against stronger tidal 
forces. In the second hand, the Gaussian nature of the fluctuation field means that a good 
approximation to the number density of bound objects that have mass m at time z is given 
by considering the barrier crossing statistics of many independent and uncorrelated random 
walks, where the barrier shape B{m,z), is connected to the collapse threshold. Moreover, 
using the shape of the modified barrier in the excursion set approach, it is possible to obtain 
a good fit to the universal halo mass function. ^ 

As previously reported, the excursion set approach allows one to calculate good ap- 
proximations to several important quantities, such as the "unconditional" and "conditional" 
mass functions. Shcth & Tormen (2002) (hereafter ST) provided formulas to calculate these 
last quantities starting from the shape of the barrier. They also showed that the "uncon- 
ditional" mass function, which is the one we need now, is in good agreement with results 
from numerical simulations. Using the barrier shape obtained in Del Popolo & Gambera 
(1998), obtained from the parametrization of the nonlinear collapse discussed in that paper, 
together with the results of ST in order to study the "unconditional" and "conditional" 
mass function. Assuming that the barrier is proportional to the threshold for the collapse, 
similarly to ST, the barrier can be expressed in the form: 



B{M) = 5, = 6o 



1 + 



rta^^_dr 



1 + 4^' 



where 5co — 1-68 is the critical threshold for a spherical model, ri is the initial radius, rta is 
the turn-around radius, L the angular momentum, ai — 0.585 and /3i = 0.46. ^ As described 
in Del Popolo (2000b), the mass function can be approximated by: 



n(m, z) ~ 1.21 



rn? dlogm 



1 + 



0.06 



,0.585 




exp {— ai/ 



1 + 



0.57 



xO.585 



/2} (2) 



where a = 0.707. Eq. (2) can be converted from a mass function to a luminosity func- 
tion. Following Mathiesen & Evrard (1998) and Reichart et al. (1999), I assume that the 
bolometric luminosities are given by: 



LboI oc MP{1 + zY 



(3) 



^Note that at present there is no good numerical test of analytic predictions for the low mass tail of the 
mass function. 

^The angular momentum appearing in Eq. (1) is the total angular momentum acquired by the proto- 
structure during evolution. In order to calculate L, I'll use the same model as described in Del Popolo & 
Gambera (1998, 1999) (more hints on the model and some of the model limits can be found in Del Popolo, 
Ercan & Gambera (2001)). 



- 7- 




43.5 44 44.5 

log Lno=i (erg/sec) 



Fig. 1.— The ROSAT BCS luminosity function. The sohd hue is the best fit of Eq. (17) 
of Reichart et al. (1999) to all 12 luminosity bins. The dotted line is the best fit of Eq. 
(13) to all 12 luminosity bins. The dashed line is the best fit of Eq. (13) to all but the 
highest-luminosity bin. 
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Fig. 2. — The 1, 2 and Sa credible regions of the posterior probabiUty distributions of the 
fit of Eq. (13) to all but the highest- luminosity bin. 
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(a) 




Fig. 3. — The marginalized posterior probability distributions P{Qm) (left panel) and P{n) 
(right panel) of the fit of Eq. (15) to the 0.14 < z < 0.6 revised EMSS clusters and the 
ROSAT BCS luminosity function. The solid line is the result of the calculation of this paper, 
the dotted fine that of Reichart et al. (1999), and the dashed fines mark the 1, 2, 3 cr credible 
intervals. 
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and that: 

Lx oc fxT-^L^oi (4) 

where fx — 0.989 ± 0.014 and /3 = 0.407 ± 0.008 for the representative temperature range 
of EMSS (3 < T < lOkeV) and fx = 1.033 ± 0.012 and P = 0.472 ± 0.008 for the represen- 
tative temperature range of BCS (1.5 < T < 12keV). In order to remove the temperature 
dependence introduced by Eq. (4), we need an M-T relation. The M-T relation that I'll 
use, obtained in Del Popolo (2000a), is based on the merging-halo formalism of Lacey & 
Cole (1993), accounting for the fact that massive clusters accrete matter quasi-continuously, 
and is an improvement of a model proposed by Voit (2000) (hereafter V2000), again to take 
account of angular momentum acquisition by protostructures. 

Under the standard assumption of the Intra-Cluster (IC) gas in hydrostatic equilibrium 
with the potential well of a spherically symmetric, virialized cluster, the IC gas tcmpcraturc- 
mass relation is easily obtained by applying the virial theorem and for a flat matter-dominated 
Universe it is given by (Evrard 1990, Evrard et al. 1996, Evrard 1997): 

T - {(^-^h'^'keV) [j^j il + z) (5) 

The assumptions of perfect hydrostatic equilibrium and virialization are in reality not com- 
pletely satisfied in the case of clusters. Clusters profile may depart from isothermality, with 
slight temperature gradients throughout the cluster (Komatsu & Seljak 2001). The X-ray 
weighted temperature can be slightly different from the mean mass weighted virial temper- 
ature. A noteworthy drawback of previous analyzes has been stressed by Voit & Donahue 
(1998) (hereafter V98) and Voit (2000) (hereafter V2000). Using the merging-halo formal- 
ism of Lacey & Cole (1993), which accounts for the fact that massive clusters accrete matter 
quasi-continuously, they showed that the M-T relation evolves, with time, more modestly 
than what expected in previous models predicting T oc {1 + z), and this evolution is even 
more modest in open universes. Moreover, recent studies have shown that the self-similarity 
in the M-T relation seems to break at some keV (Nevalainen, Markevitch & Forman (here- 
after NMF); Xu, Jin & Wu 2001). By means of ASCA data, using a small sample of 9 
clusters (6 at 4 keV and 3 at ~ 1 keV), NMF has shown that Mtot oc T^'^^^"-^^ for the 
whole sample, and Mtot oc excluding the low-temperature clusters. Xu, Jin & Wu 

(2001) has found Mtot oc T^-60±o-04 (^sing the (3 model), and Mtot oc T^-8i±o i4 ^y means of 
the Navarro, Frenk & White (1995) profile. Finoguenov, Reiprich & Boeringer (2001), have 
investigated the M-T relation in the low-mass end finding that M oc T""^, and M oc T'~3/2 
the high mass end. This behavior has been attributed to the effect of the formation redshift 
(Finoguenov, Reiprich & Boeringer 2001) (but see Mathiesen 2001 for a different point of 
view), or to coohng processes (Muanwong et al. 2001) and heating (Bialek, Evrard & Mohr 
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2000). Afshordi & Cen (2001) (hereafter AC) have shown that non-sphericity introduces an 
asymmetric, mass dependent, scatter for the M-T relation altering its slope at the low mass 
end (T ~ 3 kcV). 

Clearly this has effects on the final shape of the temperature function. In the following, I'll 
use a modified version of the M-T relation obtained improving V98, V2000, to take account 
of tidal interaction between clusters. This M-T relation is given by: 



kT ~ 8keV 



Ml 



1 I (tn\l I K{mi,x) 
mi^ \ t J M8/3 



2 

3 Ko(mi,x) 



mi \ *o / 

(see Del Popolo 2000a for a derivation), where Mq ~ 5 x 10^'^h~^MQ, ta ~ 
mi — 5/{n + 3) (being n the spectral index), and: 

K{mi,x) — Fx {mi — 1) LerchPhi{x,l,3mi/5 + 1) — 
F (mi — 1) LerchPhi{x, 1, 3mi/5) 



(6) 



ifo(l-Qi„-nA)2 ' 



(7) 



where F is defined in the Del Popolo (2000a, b) and the LerchPhi function is defined as 
follows: 



LerchPhi{z, a,v) — 



n=0 



(8) 



where x = 1 + (^)^/'^ which is connected to mass by M = A{qx ^"^I^ (V2000), and where 
K^ijni.x^ indicates that K{rni,x^ must be calculated assuming t = to- 

Eq. (6) accounts for the fact that massive clusters accrete matter quasi-continuously, 
and takes account of tidal interaction between clusters. The obtained M-T relation is no 
longer self-similar, a break in the low mass end (T ~ 3 — 4keV) of the M-T relation is present. 
The behavior of the M-T relation is as usual, M oc at the high mass end, and M (xT'^, 
with a value of 7 > 3/2 in dependence of the chosen cosmology. Larger values of 7 are 
related to open cosmologies, while ACDM cosmologies give results of the slope intermediate 
between the flat case and the open case. 

By means of the previous M-T relation, Eq. (4) reads: 



ml~^ \ t > ~^ M8/3 



J_ _|_ I ki\ ^ _|_ Koi'^-L'^) 



[i + zY 



(9) 



where p = 1.77 and s = 3.14-0.65^^, for EMSS, and p = 1.86 for BCS (Reichart, Castander 
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& Nichol 1998; Reichart et al. 1999). ^ For tn » t, Eq. (9) reduces to: 



(10) 



(see Reichart et al. 1999, Eq. 11), since the M-T relation reduces to the standard M-T 
relation (T oc M'^/^{l+z)). Taking account only the change in the mass function, substituting 
Eq. (10) into Eq. (2), defining Lx = x{z)Li (sec Reichart et al. 1999 for a definition of 
x{z)), and assuming a scale- free mass density fluctuation power spectrum of power law index 
n, so that the present variance is: 



/ , / ^ , -(3+n) 

^o(M)=a8(— ) — 



where Mg 



2G 



(8h-iMpc)3 = 6.0 X W^VL^\i-^Mq (Voit 2000) (~ W^Mq, in a flat 



universe with h — 0.5). I obtain the following luminosity function: 



^ = C2F{z, Li) exp [-G{z, Li)] 



(12) 



A(.) = c,^/^^5.o(.)(l + .)''^^ .,.-V2-^ 



/x 



: , X 1 /2 

Zp-213 x[z) ' 3p-2/: 



f2{z) = Cl 



_ -| /Q 3 — n + 6 Q-|-2 g n 



f6c. + 2an ^ . o„\, xT 1 ( - s + (3-n+6 Q + 2 g n) ] , /p 3-n+6 

3p+2f 5co(;2)^ ^ ^ (1 + ^)^" ^ J /x ^ 

. _l/2 3-n+6g+2an 

aia;(2;j ' 3p-2/3 

^/ r N r / N r 3+n (g-l)(3+u) (-1+2q)(3+u) 

G(2;,Li)= gi{z) Li3p-2 + g2{z) Li ^v-^^ ^ g^[z) LI 3p-2/3 



+6 Q+2 Q n 



-2/3 a" 



3+n 



1 3+n 

^l(z) = -Cl-3P+2f /x 3p-2/3(l+2) 



(5-/3) (3+n) ^ , ,9 , , 3+n 

3P-2/3 a5co('2) a;(z)3p-2/: 



^2(^) 

y3(^) = 



— 3 — n+3 a-|-CK n /-. \ 

= Cl -3p+2i> 



(s-/3)(3+n)(a-l) (g-l) (3+n) 

a2a;(2;) 3p-2g 



^ / N CO O ^> „ ("-!) (3+n) , , (s-/3)(3+n)( 

4o(-2)^'-'"^ /x ^-^^ (1 + ;^) ^ 

1 (3+n)(-l+2 a) /, „ x ^ / x Co /I „^ „ (-1+2 g) (3+n) . (s-ff) (3+n) (-1+2 g) 

_ci =3^+2^ a^i-^'^) <5co(-2)^'-'°^ /x ^-^^ (1+-^) — 



3p-2/ 



(-1 + 2 g) (3+n) 
02 ^(Z) 3P-2/3 



^Sincc the BCS is a local (zcff — 0.1) catalog, the value of s is unimportant. 

''The EMSS and BCS provide luminosities in their respective X-ray bands that have been computed for 
Hq = 50kmMpc~^s~^ and Om = 1- These luminosities are denoted by Li 

^as is the amplitude of the mass density fluctuation power spectrum over spheres of radius 8/i~^Mpc, 
and Ms is the mean mass within these spheres 
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where oi = 0.06 and 02 = 0.57 and the constants ci and ci corresponds to the constants a 
and c in Reichart et al. (1999), respectively. As in Reichart et al. (1999), C\ and Ci depend 
upon (jg and the factor of proportionahty of Eq. (9). 

Taking also account of the changes in the M-T relation, and using Eq. (9), the XLF 
function can be written as: 

dnc(m, z) ^ A{m,z) 

dLx B{m,z) ^ ^ 

mi' — m mi' — m mi 

+ m('-«) F ■ LPhmi) ^^^'^(-m^^^'+^+^'^+g+aa) + ^{-V+^+m^+q) _ ^{-p+l+m^+q+a^) 

+ m^-P+^+^'^+D timi - m^-p+^+m2+2a^,) p . ipj^^^ ^ 2 ^(-p+i+^a+as) p . LPhmi^ 

- F ■ LPhmi^ + p . LPhmi - 2 m(-f+^+"^2+"3) F ■ LPhmi 
+ F ■ LPhmi - m^-P+'^+"'^+^"'^"+''+'''^ ai a^"") S^'^""^ 

+ m^-P+i+m2+2m4a+q) ^(-a) ^(-2 a) _ ^{-p+l+m2+2 a+q+aa) ^(-a) ^(-2 a) 
+ ^(-P+l+"»2+2m4a+g) ^(-a) ^(-2 a) 

- ^(-P+l+"^2+2m4a+2a3) ^(-a) ^(-2 a) ^ . pp^j^^ 

+ 2m(-^'+^+"*^+^"'^"+"3) Oi a(-") 5^-2") F • LPhmi^ 

- ^(-P+l+^2+2m4a) ^(-a) 5(-2a) ^ . pPhmi^ 
+ ^(-P+l+"»2+2m4a+2a3) ^(-a) ^(-2 a) p . ^p/^ 

+ ^(-P+l+"»2+2m4a) ^(-a) ^(-2 a) p . pPhmi) 



B{m,z) — m^^'^cify^{(5Fa^mi' — (5Fa^mi — I3m''r + m'^p — pF LPhmi' + pF LPhmi + 

m'^pti mi - m''+"3p - /SFa^ mi'mn LPh - [3rFLPh)mi + Prmi'FLPh + jSFLPhqmi + 
PPas mil LPhmi + l3m"-'-^Fa^ mi + Pm"-'-^ FLPha^ mi - Pm'^rti mi - 2 m^^pF LPhmi + 
2 m^^pFLPhmi^ + 2 (3m'''FLPhmi^q - pm""' FLPhmi^as - 2 pm'^'FLPhqmi + 
2 Pm^^rFLPhmi + 2 Pm^^^Fa^ mi^mn LPh - 2 Pm^^Fa^ mn LPhmi - f^m^^Fa^ mi^ - 
f3FLPhmi^q - 2 l3m"-^rmi^FLPh + /3m*+"3^ + m^^^pFLPhmi + (im^^'^Fa^ mu LPhmi + 
(3m'^''^FLPhmi^a^ - m«+>ii mi + Pm^"' FLPhqmi + Pm^^'rmi'^FLPh - 
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(3m^''^FLPha^ rui - (3m'^''^Fa^ mi^mu LPh) 

where, in order to simplify the previous expression, I have defined: m — mi = 
5/ (n + 3), 1712 = n/G, rriz — 1 + n/3, = —1/2 — n/6, a, ai, and 02 are the same of Eq. 
(12 ), as = -5/(3ml), rriu = 3ml/5 + 1, h ^ {ta/tf^\ r = 2/3, q = 8/3, 4o(^) = S, and 
LerchPhiim"^ , l^mu) = LPh. 

Then the luminosity function can be written in implicit form as: 



dnc{m,z) _ A{m,z) , , 

dLx ~ B{m, z) ^ ' 



Lx = x{z)Li = ci/x (m''(l + mi + m-''+''^ F LPhrui^ - m~'^FLPhmi^ - m-''+''^ FLPhrui + 
m-'iFLPhmi))-'^wPml {l + z)' 

The total number of X-ray clusters observed between luminosity and redshift limits 
L\ < LI < Lu and z\< z < Zu, i.e. the cumulative luminosity function is given by: 

Ar(Li, L„ ^1, z^) = J^^ £ yl(Li, z) ^^^^^^j' dL,dV{z) (15) 

where A[Li, z) is the area of the sky that an X-ray survey samples at redshift 2; as a function 
of luminosity Li. In the case of the EMSS, this quantity is given by (Avni & Bahcall 1980; 
Henry et al. 1992, Nichol et al. 1997) 

A{Li,z)=A{F^^ = F{Li,z)) (16) 

where A{Fii^ is the area that the EMSS surveyed below sensitivity hmit F^^ (see Henry et 
al. 1992): 

k{z) ATcdl{z) 

being k{z) the -correction from the observer frame to the source frame for a T = 6keV X-ray 
cluster. A(Li,z) can be calculated similarly to Reichart et al. (1999) (see their Fig. 1). 
Finally 

. _ ^c'dz (QruZ + (n^n - 2)((Q^z + 1)^ - 1))2 
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is the comoving volume element. The previous model, has reported by Reichart et al. (1999) 
consists of the parameters: Hq, fx, P, p, s, Ci, C2, n and fl. Fixing the value of Hq, using the 
quoted values of fx, P, p, s, and noticing that the normalization parameter, C2, drops out of 
the Baycsian inference analysis of this paper (see Eq. (20)), we are left with the parameters 
fim, n, c\. It is possible to obtain credible intervals for the values of the quoted parameters 
as described in Reichart et al. (1999). As in the quoted paper, the parameter n and ci are 
constrained with the local 2;efr ~ 0.1 luminosity function of the BCS. Althought such a sample 
may not have enough redshift leverage to adequately probe Jim, its large size makes it an 
excellent sample to constraint the parameters n and c\. These constraints can be combined 
with the EMSS results to better constrain fim- So letting z — Zeg, A(Li,Zeff) = 1 (see 
Reichart et al. 1999 for a discussion), using the binned BCS luminosity function of Fig. 1 of 
Ebeling et al. (1997), I find the parameters by using a likelihood function given by e~^^/'^, 
with -3 < n < and < Ci < 10. Fitting Eq. (14) to aU data, 1 obtain n = -O.Gbtols 
and Ci = O.77I0 11- When fitting Eq. (14) to all but the highest luminosity bin, 1 find that 
n = — 1.53^o;^2 ci = l.OS^oJg. This double fitting is due to safeguard against a bias on 
n that can be introduced by the assumption z — Zes (see Reichart et al. 1999 for a deeper 
description) . 

In Fig. 1, I plot the ROSAT BCS luminosity function. The solid line is the best fit of 
Eq. (17) of Reichart et al. (1999) to all 12 luminosity bins. The dotted line is the best fit 
of Eq. (13) to all 12 luminosity bins. The dashed line is the best fit of Eq. (13) to all but 
the highest-luminosity bin. In Fig. 2, I plot the 1, 2, 3a credible regions in the n-ci plane 
for the second fit (all but the highest luminosity bin). 

As previously described, the data coming from EMSS are used to obtain constraints on 
f^m, n, Ci, by means of Bayesian inference analysis. In order to perform the quoted analysis, 
we need the posterior probability distribution for n, Ci, namely P{D,ja.,n,Ci) which is 
obtained by normalizing the product of the prior probability distribution and the likelihood 
function (Gregory & Loredo 1992). The prior assumed is fiat between < fl^n < 1-5, 
—3 < n < and < Ci < 3. In order to constraint, l^m, n and Ci, I use the likelihood 
function L{fl^,n,Ci) (see Cash 1979, Reichart et al. 1999): 



where P{Li,i , Zi\flrn, n, Ci) is the probability that the ith X-ray cluster fits our model, given 
values of fl^, n, ci. For our model, this probability is given by (e.g.. Cash 1979): 



Ntot 

L{Qra, n, ci) = Y\ P{Li,i , ZilQrn, n, Ci) 



(19) 



i=l 



P{Li,z\nrn,n, Ci) 




(20) 



N{Li, Lu, zi, z,^ ) 
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In Eq. (19), A^tot is the total number of X-ray clusters in the same region of the Li — z plane 
as that over which N{Li, L^, z\, z,u ) is defined. The value of L\ < Li < L^^ = 10^^'^erg/s. Li 
is set by the limiting flux of the EMSS: F{Li, z) = 1.33 x IQ-^^erg/cm^s (Henry et al. 1992). 
The redshift is in the range: 0.14 — zi < z < z^ = 0.6. 

The posterior probability distribution for one of the parameters, e.g., P{Qra), is ob- 
tained by marginalizing the posterior probability distribution for all the three parameters, 
P(r2jn, n, Ci), over the other parameters (Gregory & Loredo 1992). The credible regions (1, 
2, 3 a) arc determined by integrating the posterior probability distribution over the most 
probable region of its parameter space until 68,3%, 95.4% and 99.73% (respectively) of this 
distribution has been integrated. In Fig. 3, 1 plot the marginalized posterior probabil- 
ity distributions P(f2m) (left panel) and P{n) (right panel) of the fit of Eq. (15) to the 
0.14 < ^ < 0.6 revised EMSS clusters and the ROSAT BCS luminosity function. The sohd 
line is the result of the calculation of this paper, the dotted hue that of Reichart et al. (1999). 

^ In Fig. 3a-b, the dashed hues represents the 1, 2, and Sa credible intervals, the dotted 
line of Fig. 3a represents P(f2m) versus Qm, while Fig. 3b represents P{n) versus n. The 
figure shows that — l-lS^ass ^ — ~^-^^-oai- The previous result shows that the 
change in the mass function and M-T relation gives rise to an increase of and n of ~ 20%. 

The lesson from the previous calculation is that taking account of non-sphericity in 
collapse and the fact that massive clusters accrete matter quasi-continuously gives rise to a 
noteworthy change in the prediction of cosmological parameters, as D,^. Here, I am princi- 
pally interested in studying the effects of these changes on the values of cosmological param- 
eters, and not in the peculiar value obtained. In this light, the large value of fl^ obtained is 
not of particular importance or meaning, but is simply a consequence of the high value of f2m 
obtained in Reichart et al. (1999) analysis. As an example, I have also estimated the value 
of fim following Borgani et al. (2001). Analyzing the ROSAT Deep Cluster Survey (RDCS) 
and using the XLF to obtain constraints on cosmological parameters, Borgani et al. (2001), 
found that = 0.351^;^^. In their study, they used the ST mass function instead of the 
usual PS and Eke et al. (1998) M-T relation. Using their method and data, but our mass 
function and M-T relation, one obtains larger values of Jim (fim ~ 0.4 ± 0.1) that differently 
from the previous analysis (Reichart et al. 1999) exclude an Einstein-dc Sitter model. Note 
that since the RDCS data are not available to date, I used data taken from the XLF in 
Rosati, Borgani & Norman, (2002). In any case, even using all RDCS data small changes 



^Namely, the likelihood function is that of the EMSS, Eq. (19). Instead of assuming a flat priori 
distribution for all three parameters, I assume a flat priori between < flm < 1-5, and use the posterior 
probability distribution of the BCS-Pscsin, ci) (see Fig. 2)-Eq. (12), and the effective redshift of the BCS, 
ZeS — 0.1, to determine the full prior probability distribution: PBCs{^m,Ti,ci). 
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are expected, and the quoted result, namely that the Einstein-de Sitter model is excluded, 
is not changed, but a larger value of is obtained. 

3. Constraints to Om, n, and as from the XTF 

As previously reported, the mass function (MF) is a critical ingredient in putting strong 
constraints on cosmological parameters (e.g., Jim)- Observationally the local mass function 
has been derived from measuring masses of individual clusters from galaxy velocity disper- 
sions or other optical properties by Bahcall and Cen (1993), Biviano et al. (1993), and 
Girardi et al. (1998). However, the estimated virial masses for individual clusters depend 
rather strongly on model assumptions. As argued by Evrard ct al. (1996) on the basis 
of hydrodynamical N-body simulations, cluster masses may be presently more accurately 
determined from a temperature measurement and a mass-temperature relation determined 
from detailed observations or numerical modeling. Thus alternatively, as a well-defined ob- 
servational quantity, the X-ray temperature function (XTF) has been measured, which can 
be converted to the MF by means of the mass-temperature relation. 

The cluster temperature function is defined as: 

N{T,z)^N{M,z)^ (21) 

While the mass function, N{M, z), gives the mass and distribution of a population of evolving 
clusters, the Jacobian ^ describes the physical properties of the single cluster. 

Comparison of the predictions of the PS theory with the SCDM and OCDM cosmolo- 
gies, performed by Tozzi & Govcrnato (1998) and Governato et al. (1999), have shown 
discrepancies between PS predictions and N-body simulations, increasing with increasing z. 

Before going on, I want to discuss about the use of Eq. (2). Jenkins et al. (2001), 
obtained a mass function, which is regarded as perhaps the most accurate model to date, 
from the Hubble volume simulations. 

Jenkins et al. (2001) showed that, although the mass functions in their simulations 

scaled in accordance with the excursion set prediction, Shcth & Tormcn 2002 (Eq. (2)) 
slightly overestimated the unconditional mass functions in their simulations. It interesting 
to note that as shown by Pierpaoli et al. (2001), using Sheth & Tormen (2002) or Jenkins et 
al. (2001, 2003) mass function, rather than PS, the value of parameters like erg are changed 
by a small amount (4-8 %). Sheth & Tormen (2002) also showed that changing the parameter 
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Fig. 4. — The 68% confidence contours for the parameters n, and for the open modeL 
The dashed hnes are hnes of constant a%. 
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Fig. 5. — A(likelihood) for the parameter The sohd hne is obtained from the model 
of this paper while the dotted line is that calculated by Henry (2000). The dashed hnes 
represent various confidence levels (65%, 90%, 95%, 99%). 
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Fig. 6. — The 68% confidence contours for the parameters Ug, and for the open model 
(see also Henry (2000), Fig. 9). 
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Fig. 7. — Constraints on and obtained using the same 25 clusters used in Henry 
(2000), for the local sample, while the high redshift sample is constituted from all the EMSS 
clusters with z > 0.3 and RX J0152. 7-1357 (see Henry 2002). The solid lines are the 1 and 
2 a contours obtained using the mass function and the M-T relation of this paper, while the 
dashed line is the 1 o" contour obtained using the PS mass function and the M-T relation of 
Pierpaoh et al. 2001. 
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Fig. 8. — (a) Comparison of various mass functions. The dotted line represents Sheth & 
Tormen (2002) prediction, the sohd hne that of Jenkins et al. (2001) and the dashed hne 
that of Del Popolo (2000b). 



-24- 



a in their Eq. (2) from 0.707 to 0.75 reduces the discrepancy between it and the simulations 
substantially. 

In Del Popolo (2002a), I showed that Eq. (2) of this paper, in agreement with Jenkins 
et al. (2001), predicts smaller values of the mass function expecially for high with respect 
with Sheth & Tormen's predictions. In other terms, Eq. (2) of this paper is in very good 
agreement with Jenkins et al. (2001). In Fig. 8 I plot a comparison of the various mass 
functions: the dotted line represents Sheth & Tormen (2002) prediction, the solid line that 
of Jenkins et al. (2001) and the dashed line that of Del Popolo (2000b). As shown, Jenkins 
et al. (2001) mass function is almost indistinguishable from that used in this paper. For 
what concerns small masses, our formula, in agreement with Sheth & Tormen 2002, differs 
dramatically from the one Jenkins et al. (2001), propose. Simulations currently available do 
not probe the regime where u < 0.3 or so (the Jeans mass is at about u < 0.03 (Sheth & 
Tormen 2002)). New simulations are needed to address which low mass behavior is correct. 
In other terms, the mass function obtained in this paper is in very good agreement with 
Jenkins et al. (2001) in the regime probed by simulations. Moreover, the constraints obtained 
in almost all the papers in literature used the PS mass function except a few papers (e.g. 
Borgani et al. 2001; Henry 2002) and the M-T relation the usual one obtained from the virial 
theorem. In other words, this paper introduces noteworthy improvements on the previous 
calculations in hterature. The mass entering in Eq. (2) should be interpreted as the mass 
contained inside a radius, ri8o encompassing a mean overdensity p — 180p. However, scahng 
relations connecting mass X-ray observable quantities may provide the mass at different 
values of p/p. In this case we follow White (2001) and rescale the masses assuming an NFW 
(Navarro, Frenk & White 1996) profile for the dark matter halo with a concentration c = 5 
appropriate for rich clusters (see also Pierpaoli et al. 2001, 2003; Schuecher et al. 2002 for 
more details). 

Similarly, it necessary to give arguments to use Eq. (6) instead of the new mass/X-ray 
temperature relations obtained from simulations or Chandra data within the last year (see, 
e.g., Pierpaoli et al. 2001, 2003 for a reference). As shown in Del Popolo (2002a), Eq. 
(6) reduces to a similar equation to that used in Pierpaoli et al. (2001) (Eq. 13), in the 
early-time, namely: 

M (X T3/2p-i/2A-V2 0, r^/2(A,E2)-V2 (22) 

where E{zY = Qui^ + z)^ + ^Aj and the term depending on in Eq. (13) of Pierpaoli et 
al. (2001) (which is a correction to the virial relation arising from the additional potential 
in the presence of A) is neglected since it produces only a small correction (see Pierpaoli 
et al. 2001). I also want to add that Eq. (13) of Pierpaoh et al. (2001) or Eq. (4) of 
Pierpaoli et al. (2003), comes from rather simplistic arguments (dimensional analysis and an 
assumption that clusters are self-similar) and is a good approximation to both observations 
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and simulations but this last are sufficiently computationally demanding that they cannot 
explore parameter space efficiently and so it is necessary to determine coefficients by means of 
simulations, while scalings are taken from simple theoretical models (Pierpaoli et al. 2001). 
Eq. (13) of Pierpaoli et al. (2001), is valid for systems hotter than about 3 kcV. We know 
that recent studies have shown that the self-similarity in the M-T relation seems to break at 
some keV (Nevalanien et al. 2000; Xu, Jin & Wu 2001): Afshordi & Cen (2001) has shown 
that non-sphericity introduces an asymmetric, mass dependent, scatter for the M-T relation 
altering its slope at the low mass end (T ~ 3 keV). These effects are taken into account 
by mine M-T relation, which as previously told gives same results of that of Pierpaoli et al. 
(2001, 2003) in the range of energy in which this is valid. 

In the following, I'll use the mass function modified as described in the previous section 
and an improved form of the M-T relation in order to calculate the mass function. 

Introducing Eq. (2) into Eq. (21) and using the M-T relation, in the peculiar case that 
the variance is given by: 

a = C'(^-)-(3+")/6 = Cm-(^+")/6 (23) 

I get: 



dN_ ^ N,{m,z) 

dT N2{m,z) ^ ' 



Ni{m,z) = l/2mi (3 + n)p5 Va\/2Am™2+« 

- -4^1 m"'3a52-ai-«52-2 "m'"3 (Cm'"4 )2 «a2-l/2 a^-^ "5'^-^ °m'"3 (Cm™4 )* "02^ 



(l + aia-"r2"(Cm"*^)'") 



N2{m, z) = m^/^(3 Frui^LPh qrrf^ - 3 Fmi Og mn LPh rrf^ - 2 FLPh mx^rrf^ + 

3 Fmi^as mu LPh rrf^ + 2 FLPh mi m"^ - 3 Fmi LPh qrrf^ - 3 Frm^a^ LPh m"^ + 
3 Fmi as LPh m"^ - 3 Fmi^LPh g - 3 Fnii a?, + 2 FLPh mi^ - 2 m« - 2 ti mi - 
3 Fmi^aa mu LPh + 3 Fmi 03 mu LPh + 3 Fmi^a^ - 2 FLPh mi + 3 Fmi LPh q)C ^/7^ 
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kT ~ SkeVrri^/^ 



1 I +3 _|_ -^'("'1,3:^) 
^1 (Mi5m)8/3 



mi 



mi + ^1 + M«/^ 



In order to use the same notation and variance of Henry (2000), the constant C is 
defined as: 



C = 0.675 , 



r(3 + n) sin(l/2 nvr) /^^^ ko'^Mpc 



2"n (2 + n) (1 - n) (3 - n) 



3 _i/2-l/6n 
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and then (Tg = (T(f2jn, M = 0.594 x lO^^/i^^Om), and the XTF depends on the parameters 
n, ko and The data that shall be fitted to the theory previously described, are those 
described in Section. 2 of Henry (2000). I use a maximum likelihood fit to the unbinned data 
in order to determine various model parameters. The method is described in Marshall et al 
(1983). The hkelihood function is given by their Eq. (2), which for the present situation is: 



N 



5= -2 J] In 



dzdVt 



+2 



/ dkT I 



n{Vtrni kT)Q{z, kT) 



dzd,^ 
(26) 

(Henry 2000), where N is the number of clusters observed, n{VL.^^z^kT) is the temperature 
function, Vt{z, kT) is the solid angle in which a cluster with temperature kT at redshift z 



dz 



could have been detected (the selection function) and 



dzdU 



is the differential volume. 



which is given in the Appendix of Henry (2000) (see also Henry 2000 for a description of 
^d?dn^^^ ^^^^ term of Eq. (26)). The best estimates for the parameter are obtained 

minimizing S and the confidence regions are obtained as in Henry (2000) by noticing that S 
is distributed as with the number of degree of freedom equal to the number of interesting 
parameters (Avni 1976). The model is accepted or rejected according to what prescribed in 
Marshall et al. (1983) (see also Henry 2000). 

At this point, we can fit the data described in Section. 2 of Henry (2000) to the 
theory previously described using the quoted maximum likelihood method. The most general 
description of the results requires the three parameters of the fit. I show these results in 
Fig. 4, where I plotted the results for the open model. It is straightforward to read off the 
value of n, which is at the 68% confidence for the open model. These values shows that 
the correction introduced by the new form of the mass function and M-T relation gives rise 
to higher values of (fi^ = 0.6 ± 0.13, while it is fi^ = 0.49 ± 0.12 for Henry (2000)) 
and n = -1.5 ± 0.32 (n = -1.72 ± 0.34 in Henry (2000)). The presentation in Fig. 4 is 
somewhat difficult to appreciate, so we also give the constraints for fewer parameters. In 
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Fig. 5, I plot A (likelihood) for the parameter Jim- The solid hne is obtained from the model 
of this paper while the dotted line is that calculated by Henry (2000). The dashed lines 
represent various confidence levels (68%, 90%, 95%, 99%). Constraints are relatively tight 
when considering this single parameter. Wc find that Q^^^ = O.GIq'}^ at the 68% confidence 
level and at Qni = 0.61q;2'^ the at 95% confidence level for the open model. 

The constraints in Fig. 4 arc plotted into a more conventional format in Fig. 6. Three 
parameters are still required, but the constraints on n and fcg are collapsed into (Xg . Fig. 6, 
plots the 68% confidence contours for the parameters (jg, and for the open model (see 
also Henry (2000), Fig. 9). In Fig. 7, 1 plot the constraints on Q/y and Q^i obtained using the 
same 25 clusters used in Henry (2000), for the local sample, while the high redshift sample is 
constituted from all the EMSS clusters with z > 0.3 and RX J0152.7-1357 (see Henry 2002). 
The sohd fines are the 1 and 2 a contours obtained using the mass function and the M-T 
relation of this paper, while the dashed line is the 1 a contour obtained using the PS mass 
function and the M-T relation of Pierpaoli et al. 2001. 

For a CDM spectrum, the expression for the XTF is much more complicated. It can 
be obtained combining Eq. (2), Eq. (21), and our M-T relation. The mass variance can be 
obtained once a spectrum, P{k), is fixed, by: 

a\M) = 7^ / dkk^P{k)W\kR) (27) 
27r Jo 

where W{kR) is a top-hat smoothing function: 

W{kR) = (sin kR - kR cos kR) (28) 

and the power spectrum P(k) — Ak"'T'^{k) is fixed giving the transfer function T{k). The 
CDM spectrum used in this paper is that of Bardeen et al. (1986) (equation (G3)), with 
transfer function: 

T{k) = ^^^Wf^^"^^^ ■ [1 + 3.89g + (16.1g)2 + {5A6qf + {6.71)^'/' (29) 

Bardeen et al. (1986)(Eq. (G3)), where q = q^^^^^-^- Here 6 = /)cr/(l-68p^,) represents 
the ratio of the energy density in relativistic particles to that in photons {9 = 1 corresponds 
to photons and three flavors of relativistic neutrinos). Using the data used by Eke et al. 
(1998) (see their Sec. 3 and the a maximum likefihood parameter estimation (see Eke et 



^°They combined the temperature data for 25 local clusters by Henry & Arnaud (1991) with the sample 
of 10 EMSS clusters at 0.3 < < 0.4 by Henry et al. (1997) 
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al. 1998, Sec. 4): 



N 




InL 



adzdT 



(30) 



i=l 



where 



a{z, t) 



n{z,T,p) 



dV{z,p) 



(31) 



dz 



where p represents the three model parameters being investigated, n{z,T,p) is the comoving 
cluster number density given by Del Popolo (2000b) model for the mass function, V is the 
comoving volume and q is the product of the fraction of the sky surveyed and the estimated 
completeness of the survey (see Eke et al. (1998), Sec. 4, for details). 

The results are plotted in Fig. 9. It plots the A(— Inlikelihood) for fl^ (left panel) and 
CTg (right panel), marginalized over the other two parameters. Solid lines are the prediction 
obtained from the model of this paper while the dotted hues those obtained from Eke et al. 
(1998). Dashed hnes are 1 and 2 a significance and the 3 a corresponds to the top of each 
panel. Fig. 9 shows that Qm (left panel) and erg (right panel) are increased with respect to 
Eke et al. (1998) prediction: while in Eke et al. (1998) = 0.52;|]:}^, and ag = O-GStom, 
I find that fl^ = 0.62+°;^^ and as = 0.73^q'^1. This shows again an increase in fl^, also in 
qualitative agreement with Eke et al. (1998) calculation taking account of changes in the 
threshold for collapse suggested by Tozzi & Governato (1998). 



In this paper, I have revisited the constraints obtained by several authors (Reichart et 
al. 1999; Eke et al. 1998; Henry 2000) on the estimated values of f^m, n and cg in the light 
of recent theoretical developments: new theoretical mass functions, a more accurate mass- 
temperature relation, also determined for arbitrary and Qa- Using the mass function 
derived in Del Popolo (2002b), the M-T relation derived in Del Popolo (2002a), and following 
Reichart et al. (1999), I calculated the XLF which can be apphed to the high-redshift X-ray 
cluster luminosity catalogs to constraints cosmological parameters, namely in this case 
and n. This luminosity function was applied, for a fixed value of Hq — 50kmMpc~^s~^, 
to broad subsets of the revised EMSS X-ray cluster subsample of Nichol et al. (1997) and 
to ROSAT BCS luminosity function of Ebeling et al. (1997) to constraint fim- For the 61 
revised EMSS clusters, with 0.14 < z < 0.6, I find that = l.lbtois and n = -l.bbto.tl- 
The previous result shows that the change in the mass function and M-T relation gives rise 
to an increase of and n of ~ 20%. Then, taking account of the fact that massive clusters 
accrete matter quasi-continuously and taking account of non-sphericity in collapse changes 



4. 



Results and discussion 
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the values of the estimated parameters. The principal interest in this paper is to study the 
effects of improvements on the mass function and M-T relation on the values of cosmological 
parameters, and not in the peculiar value obtained. So, I have also estimated the value 
of fl^ following Borgani et al. (2001). In their study, they used as mass function the ST 
mass function instead of the usual PS, (differently from Reichart et al. (1999)), which is a 
good approximation to the function I used in this paper. They used Eke et al.(1998) M-T 
relation, which is different from the one I used.. Using their method and data (RDCS), but 
our mass function and M-T relation, I obtained a larger value of Qm (f^m = 0.4 ±0.1 instead 
of = 0.351q;iq) that differently from the previous analysis (Reichart et al. 1999) exclude 
an Einstein-de Sitter model, but shows again that a change in the mass function and M-T 
relation influences the value of the parameters to constraint. Another possibility to constraint 
cosmological parameters is to use the XTF. 1 have repeated the Henry (2000) analysis but 
differently from the quoted paper, 1 changed the mass function and M-T relation, adopting 
again those of Del Popolo (2000a,b). The qualitative result is similar to the previous one, 
using the XLF, namely the values of the obtained cosmological parameters are modified 
by a ~ 20%. The new form of the mass function and M-T relation gives rise to higher 
values of {Q^ — 0.6 ± 0.13 in my estimation, while it is Jim = 0.49 ± 0.12 for Henry 
(2000)) and n = -1.5 ± 0.32, in my estimation, while n = -1.72 ± 0.34 in Henry (2000). 
I have also obtained some constraints on Q\ and flrn obtained using the same 25 clusters 
used in Henry (2000), for the local sample, while the high redshift sample is constituted 
from all the EMSS clusters with z > 0.3 and RX J0152. 7-1357 (see Henry 2002). The 
1 a contours obtained using the mass function and the M-T relation of this paper, and 
plotted in Fig. 7, shows that for A = 0, it is 0.32 < fl-^ < 0.57 in the case of Henry 
(2002) and 0.43 < < 0.73 in my estimation. The figure shows the constraints to 
for different values of Qa. Similar results to that obtained in the previous comparison with 
the Henry (2000) results are obtained changing data and method, by following Eke et al. 
(1998). I obtain a value of = 0.62+°;}^, while in Eke et al. (1998) fl^ = 0.52+°g^^ and 
as = 0.73lo!o6 while in Eke et al. (1998) as = 0.63^0 05. This shows again an increase in 
Qui! also in agreement with Eke et al. (1998) calculation taking account of changes in the 
threshold for collapse suggested by Tozzi & Governato (1998). As previously told, this 
paper has its aim that of studying how "systematic uncertainties" can influence the values 
of some cosmological parameters. It is well known that in literature the values obtained for 
Qra span the range 0.2 < Qm < 1 (Reichart et al. 1999). 

Sadat, Blanchard & Oukbir (1998) and Reichart et al. (1999) Blanchard & Bartlett 



^^Replacing 6c with the deS of Tozzi & Governato (1998), produces quahtatively similar changes on the 
mass function as those obtained using the Del Popolo (2002b) mass function. 
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(1998) found results consistent with — 1. Viana & Liddle (1999) found that Qm = 0.75 
with > 0.3 at the 90% confidence level and ~ 1 still viable. Blanchard, Bartlett 
& Sadat (1998) found almost identical results {Q^ ~ 0.74 with 0.3 < < 1.2 at the 
95% confidence level). Eke et al. (1998) found = 0.45 ± 0.2. It is interesting to note 
(as previously mentioned) that Viana & Liddle (1999) used the same data set as Eke et 
al. (1998) and showed that uncertainties both in fitting local data and in the theoretical 
modeling could significantly change the final results: they found — 0.75 as a preferred 
value with a critical density model acceptable at < 90% c.l. 

Different results were obtained by Bahcall, Fan & Cen (1997) (r2m = 0.3 ± 0.1), Fan, 
BahcaU & Cen (1997) (fi^ = 0.3 ± 0.1), BahcaU k Fan (1998) {n^ = 0.2+^;^) and several 
other authors. 

The reasons leading to the quoted discrepancies has been studied in several papers (Eke 
et al. 1998; Reichart et al. 1999; Donahue & Voit 1999; Borgani et al. 2001). According 
to Reichart (1999) unknown systematic effects may be plaguing great part of the quoted 
results. A list of these last is reported in the introduction to this paper. Our analysis 
shows that improvements in the mass function and M-T relation increases the value of ft^. 
The effect of this increase is unable to enhance significantly the probabihty that fl^ — 1 
in the case of constraints like that of Fan, Bahcall & Cen (1997) {Q^a = 0.3) or Bahcall & 
Fan (1998) = 0.2), and can give a small contribution even in the case of larger values 
for the value of the constraints obtained. However, in any case it shows that even small 
correction in the physics of the collapse can induce noteworthy effects on the constraints 
obtained. Moreover, even changing the data or the way they are analyzed gives different 
results. As an example, changing their low-redshift sample, Donahue & Voit (1999) showed 
a change in Q,^ from 0.45 to 0.3. Furthermore, as observations are reaching the first epoch 
of cluster assembly, treating them as dynamical relaxed and virialized systems is undoubtly 
an oversemplification. Hierarchical clustering scenario predicts that a fraction between 0.3 
and 0.6 of the z = 1 population of clusters are observed less than 1 Gyr after the last major 
merger event and then are likely to be in a state of non-equilibrium. Although the quoted 
uncertainties has been so far of minor importance with respect to the paucity of observational 
data, a breakthrough is needed in the quality of the theoretical framework if high-redshift 
clusters are to take part in the high-precision-era of observational cosmology. 
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